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The recently discovered charmed-strange meson D*q (2317) has been speculated to be a tetraquark 
mesonium. We study this suggestion with overlap fermions on 2+1 flavor domain wall fermion 
configurations. We use 4-quark interpolating operators with Z4 grid sources on two lattices (16^ x 
32 and 24^ x 64) to study the volume dependence of the states in an attempt to discern the nature 
of the states in the four-quark correlator to see if they are all two-meson scattering states or if one 
is a tetraquark mesonium. We also use the hybrid boundary condition method for this purpose 
which is designed to lift the two-meson states in energy while leaving the tetraquark mesonium 
unchanged. We find that the volume method is not effective in the present case due to the fact that 
the scattering states spectrum is closely packed for such heavy states so that one cannot separate 
out individual scattering states since the volume dependence is skewed as a result. However, the 
hybrid boundary condition method works and we found that the four-quark correlators can be 
fitted with a tower of two-meson scattering states. We conclude that we do not see a tetraquark 
mesonium in the D*q(2317) meson region. 
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1. Introduction 

A charmed-strange meson D*q(2317) has been found in recent years ^. Its mass is 
2317.8 ±0.6 MeV, which is significantly lower than the one predicted from the quark model 
This discrepancy has prompted speculations that D*q(2317) is a DK molecule [^], a four quark 
state or a dynamically generated bound C5 state through the coupling to the nearby DK thresh- 
old 

Lattice QCD is a powerful tool used to calculate the hadron spectrum from first principles. 
Some earlier works were done with dynamical lattice QCD for the meson |j8|, ^ [T^, |lT]], but the 
puzzle is not settled since the results are significantly different than the experimental mass. A recent 



calculation [12] with overlap fermion on 2 + 1 flavor domain-wall sea gives a mass of 2304 it 22 
MeV which is consistent with the experimental value. Presumably, this is largely due to the fact 
that the overlap fermion has much smaller 0{m^a^) error than the other fermion formulations and 



is suitable for both light and charm quarks []13[]. This result suggests that the D*q(2317) is just the 
scalar cs meson. 

In this article, we shall study this meson to see if it can be a tetraquark mesonium as suggested 
in the literature [||, ^. We calculate the spectrum with 4-quark interpolating operators and look for 
evidence of a 4-quark state. 

2. Strategies for the computation 

2.1 Overlap valence fermion on domain-wall sea 

The overlap fermion action obeys chiral symmetry at finite lattice spacing and is, thus, free 
of 0{a) errors. It is shown that the effective quark propagator of the massive overlap fermion has 
the same form as that of the continuum [|T^. The 0(m^a^) error, which is important in the charm 



region, is estimated to be small on quench lattices [15, 13 1 and even smaller on the dynamical 



domain-wall sea [ |12| ] due to HYP smearing. Thus, it is shown that it can be used for both charm 
and light quarks on the 24^* x 64 lattice DWF configurations [16]. 



We use the valence overlap fermion on A'y- = 2 + 1 domain wall dynamical configurations from 
the RBC and UKQCD collaborations [[l7|]. This is a mixed action approach. Since the valence is 



a chiral fermion, only one extra low-energy constant Amk needs to be determined [18, 19], which 



turns out to be small [16]. 



2.2 The Z4 grid source 

We introduce Z4 grid sources to gain more efficiency. A Z4 grid source is defined as: 

T](x) = ^e(r)% (2.1) 

where is a sparse grid of lattice sites on timesUce t = 0, and (/) G { 1 , /, — 1 , — /} is the Z4 random 
phase on site /. 



* Speaker. 

1"For the XQCD Collaboration. 



2 



Z)*q(2317) with chiral fermions 



M. Gong 





Figure 1: A sketch of the correlators with Z4 grid sources, a) Left panel represents the first term in Eq. ( 2 .4 ) 



which is a signal, b) Right panel represents the second term in Eq. (2.4) which is a noise. 



The corresponding quark propagator is: 

G{y,^) = D-\y,x)ri{x) 

= '£d{r}D-\y,x)5,j='£d(i)G{y,r) 



(2.2) 



And the anti-quark propagator is: 

G{7l,y) = {D^-\y,x)7l*{x))*'''-^' 

= £0(7)G(7j) 



(2.3) 



where the Td and Tc are transpose operations on dirac space and color space respectively. 

On space-time dimensions with implicit color and dirac indices, all connected 4-quark corre- 
lation functions have the form: 

C(y,77) = rr(riGi(j,77)r2G2(j,77)r3G3(77 J)r4G4(T] J)) 

= £ d{h)e {i2 )e{h)e {h)Tr (y,Gi {y, h )r2G2 {y, r2)r3G3 (fa , 3;)r4G4 ih , y) ) 

= Tr(riGi{y,T)r2G2iy,l)r3G3(i,y)r4G4(i,y)'^ 

+ £ Q{ii)Q{i2)Q {i'3 ) e (r4 )Tr(r,Gi (y, h )r2G2 (y, r2)r3G3 (73 , j)r4G4 (74 , y) ) 

(2.4) 



where the F's are 7 matrices. 

The second term is stochastically eliminated with many gauge configurations and/or many 
noise sources. In this limit, the correlator is a good approximation of the sum of many correlators 
with point sources. Figure ^ shows the different diagrams of the two terms. 

2.3 Hybrid spatial boundary condition 

Both 2-meson states and the possible 4-quark state are present in Eq. ( |2.4| ) and it is difficult to 



distinguish them, so we adopt the hybrid spatial boundary condition to help in this regard [21 1 
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If we use the periodic spatial boundary condition, the available momenta on quark propagators 
arem the set P = {0 , ± ^ , ± ^ , . . . } . On the other hand, if we use the anti-periodic spatial boundary 
condition, the available momenta on quark propagators are in the set A = {ib£,ib^,...}. We 
calculate the quark propagators with the periodic spatial boundary condition while calculate the 
anti-quark propagators with the anti-periodic spatial boundary condition. With this hybrid spatial 
boundary condition, the momenta of mesons are each in the set P-l-A = A while the momenta of 
4-quark state are in the set 2P + 2A = P. The lowest momentum of a 4-quark state is (0, 0, 0), while 
the lowest momentum of a qq meson state is lifted to (±^, it^, ib£). 

The D;o(2317) state is about 45 MeV lower than the DK threshold. If D;()(2317) is a tetraquark 
mesonium, it would remain at its position when the hybrid boundary condition is imposed as com- 
pared to the periodic condition. On the other hand, the DK threshold will be shifted up by 177 
MeV on the 24^ x 64 lattice, making it easier to identify the tetraquark mesonium state. 

3. Simulation results 

3.1 Simulation parameters 

We use the Nf = 2 + \ domain wall dynamical configurations from RBC and UKQCD collab- 
orations. The lattice sizes are 16^ x 32 and 24^ x 64 with lattice spacing 1/a = 1.73(3) GeV. For 
each ensemble, we use 37 configurations. 

The propagators are calculated with the multimass algorithm with masses ranging from 0.00140 
to 0.70. For this case, the mass oiujd quarks is set to 0.0135, the corresponding is 310 MeV 
which is near the ujd quark mass in the sea. The strange quark mass is set to 0.067 and the charm 
quark mass is set to 0.68. The grid spacing of the source is set to 8 on spatial dimensions. We have 
one Z4 grid source per configuration. 

3.2 Tlie interpolating operators and possible states 

The 4-quark interpolating operator we use is c^'yiu^u^y^s^ . It has the quantum number of 
D*o(2317) except with isospin 7=1. 

The correlation functions contain DK scattering states, DgK scattering states, a possible tetraquark 
state, and other higher states. 

The data are not good enough to determine the ratio of the spectral weights of the DK scattering 
states and the D^n scattering states exactly and the fitting results suggest that the spectral weights 
of DK states and the D^n states are of the same order. Therefore, we assume that all the scattering 
states have the same spectral weight. We propose to check the assumption carefully with open-jaw 
diagram calculations in the future. 

The length of the time dimension of the 16^ x 32 lattice is relatively small, and therefore the 
wrap-around states cannot be neglected. If the meson-meson interaction is neglected, the correla- 
tion functions of the wrap-around states have the form: 

C(/) = ^N ^^-E^t-E^ij-t)^ ^^ y^ 

In this case, the spectral weights of the forward propagating 2-meson scattering states and the 
wrap-around states are the same. Considering the states and the assumptions mentioned above, we 
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Table 1: Fitted results 



Lattice 


xVd.o.f. 




W 


E'a 


W 


16^ X 32 Periodic 


0.463 


0.0040(64) 


1.07(6) X 10-^ 


2.26(3) 


1.83(24) X lO"'* 


16^ X 32 Hybrid 


1.074 


0.0037(64) 


2.02(1) X 10-' 


2.23(9) 


1.71(82) X 10-"^ 


24^ X 64 Periodic 


1.264 


-0.0019(40) 


1.74(1) X lO"'' 


2.18(3) 


2.83(32) X 10-4 


24^ X 64 Hybrid 


0.946 


-0.0072(41) 


1.87(9) X 10-' 


2.27(9) 


3.62(166) X 10-4 



constmct a model for the fitting of our data: 

C(0 =Cz),.(0+Cz)^(0+C^J(0+C^r(0+Cr(0 (3.2) 

In this model, the 2-meson scattering states, the 2-meson wrap-around states, and the 4-quark state 
are: 

-[EA{p)+EB{-p)+^E)t 

Cab = W £ —^-^-^^ + (t^T-t) (3.3) 
pep/A Ea{p)Eb{-p) 

-EA(p)t-EB(-p)(T-t) 

pep /A Ea{p)Eb{-p) 

Crit) = W'e-^'' (3.5) 



where the "AB" can be "DK" or "D^n" and the particle energy of A or B is Ea/b{p) = y "^a/b + P^- 
The momentum p is summed over the corresponding set of momenta. We have calculated with both 
the original periodic spatial boundary condition and with the hybrid spatial boundary condition. 
In the model, the qq masses are input from fits of the single meson correlation functions. We 



have four parameters to fit in the model. The W in Eqs. (3.3) and (3.4) is the spectral weight of 



all the scattering states and wrap-around states, the AE in Eq. ( [3.3| ) is the mean effective energy 



of meson-meson interactions, and E' and W in Eq. are the mass and the spectral weight of a 
possible tetraquark state respectively. 



3.3 Fitted results 

The fitted results are tabulated in Table HI 

We note first that AE is consistent with zero, which suggests the meson-meson interactions 
do not play an important role and can be neglected. We note that on the 24^ x 64 lattices, the W 
fitted from hybrid boundary condition coincides with that from periodic boundary condition. This 
confirms that our assumption on the fitting model are reasonable. On the 16^ x 32 lattices, the W 
fitted from periodic and hybrid boundary condition do not agree with each other. We suspects that 
it is due to the finite volume effect with smaller time extent. 

The fitted E' are all larger than 3.7 GeV. We tried to fit it with a shallow constraint around 
2.3 ^ 2.4 GeV, but we cannot find a state there. Therefore, the fitted E' shows the effects of higher 
states and there is no 4-quark state near £)*q(2317). 
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Ratio of the correlation functions (0,33/0243) 



Ratio of the correlation functions (C^g3/C243} 









y D-K states only 




/ Ds-jt states only 




/ D-K and Ds-ji states 




/ 4-quark state only . 




' Simulation data 1 — 1 — 1 















D-K states only 

Ds-7t states oniy 

D-K and Ds-ti states 

4-quark state oniy . 

. Simuiation data 1 — 1 — 1 



Figure 2: The ratio of the correlation functions (Cigs / C243) with periodic/hybrid boundary condition. 



3.4 The volume dependence 



Another method [ |22| ] to distinguish the tetiaquaik mesonium and the scattering states is the 
volume dependence. The spectral weight of a single-particle state has little volume dependence 
while that for a two-particle scattering state exhibits a behavior. 

We have tried this method and find that it is not a good criterion for our case. In the charm 
region, there are towers of scattering states with small energy gaps. We cannot separate the states 
and the volume dependence of the correlation functions becomes very complex. Figure ^ shows 
the volume dependent behaviors of many states. The horizontal axis is time and the vertical axis is 
the ratio of the correlation functions on different volumes. The lines show the behaviors of a tower 
of DK states, a tower of D^tt states, the sum of DK states and DsTt states, and a tetraquark state 
only. The points with error-bars on the plots show our real data. 

For periodic boundary condition case, the tetraquark state should have the ratio 1 on all time 
slices while a scattering state should have a ratio at 3.4. If there are many states in the correlation 
functions, the ratios run from 1 to 3.4 while the time increases. It shows that we need a long time 
extent to observe the \/Lr' behavior of scattering states. For the hybrid boundary condition case, 
the conclusion is similar. 

The real data, however, contains DK states, D^tt states, other scattering states and excited 
states. Therefore, the behavior is more complex and it is difficult to find a simple model to fit it. 
We conclude that the time extent of our lattice is not long enough to resolve individual scattering 
state. This renders useless the volume study to discern the number of particles of the state. 

It is claimed in the lattice calculation [20] that tetraquark mesoniums are found on (1.8)^ x 3.5 



fm and (2.1)^ x 4.2 fm lattices through the volume study. The time extent is about the same as that 



of our lattice. In view of this study, the time extent in Ref. [gOp is similarly shorter than needed 
to resolve the tower the scattering states and, thus, suffers the same shortcoming as revealed in the 
present study. 



4. Conclusions 



We have studied the scalar charmed-strange meson with chiral fermions. The correlation func- 
tions with the 4-quark interpolating operator are fitted and no 4-quark state is found in our data. 
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The result shows that the meson D*q(2317) is not hkely a 4-quark state. 

The fact that D*q(2317) mass is well reproduced as a cs meson in the study with chiral 



fermions [|12|], the present calculation reinforces this interpretation. 

We find that the volume method is not effective in the present case due to the fact that the 
scattering states spectrum are closely packed for such heavy states and one cannot separate out 
individual scattering states and its volume dependence is skewed as a result. 

The open-jaw diagram will be studied in the future. More information about the spectral 
weights can be extracted in the study and we can thus make a more precise fitting model. 

The 7 = channel will be studied with disconnected diagrams and variational method. 

With similar methods, other interesting states in the charm region can be studied, such as 
X(3872),Z+(4430), etc.. 
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